Blind analysis results of the TWIST experiment 



Anthony Hillairet, for the TWIST Collaboratior{3 



University of Victoria, Victoria, Canada 

Abstract. The TRIUMF Weak Interaction Symmetry Test (TWIST) experiment was 
designed to test the standard model at high precision in the purely leptonic decay of polarized 
muons. A general four-fermion interaction model is used to describe the muon decay. TWIST 
measures three of the four muon decay parameters of this model, p, S and P^^, from the shape of 
the momentum-angle spectrum. The results of this model independent approach are compared 
to the standard model predictions and used to constrain new physics. 

Our collaboration has finalized the blind analysis of the final experimental data taken in 
2006 and 2007. This analysis mostly reached our goal of a precision of an order of magnitude 
improvement over the pre- TWIST measurements. 



1. Introduction 

The discovery of new physics is expected at the high energies reached by the LHC. However 
low energy physics such as muon decay can also play a complementary role in probing physics 
beyond the standard model (SM). 

Muon decay is a purely leptonic process that is ideal for testing the weak interaction at 
high precision. A model independent approach is possible due to the large mass of the W 
boson compared to the muon mass. The most general Lorentz-invariant, derivative-free, lepton- 
number-conserving matrix element M describing muon decay can be written in terms of helicity- 
preserving amplitudes as [1] 



where g^j are the complex weak coupling constants and are the possible interactions (scalar, 
vector, tensor). In this notation, the SM postulates that gYi = 1, and g'^j = otherwise. If 
the polarization of the decay positron is undetected, then the differential decay rate can be 
expressed as 
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i=L,R 
j=L,R 
K=S,V,T 
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Figure 1. The TWIST spectrome- 
ter. A highly polarized muon beam 
is stopped in a thin metal target. 
The muons decay to positrons that 
are tracked in a 2 T magnetic field 
by an array of planar drift cham- 



bers. 



where x = Ef./E^a_^,, -Emax. is the maximum energy of the positron, 6 is the angle between the 
muon polarization and the positron momentum, = \P^\ (the degree of muon polarization), 
and 
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The R.C. terms are radiative corrections, which become more significant as x approaches one. 
The dimensionless electron mass is defined by xq = m-e/i^max. • The muon decay parameters /?, 
5, ^ and r] are bilinear combinations of the weak coupling constants gij. The TWIST experiment 
measures p, 5 and ^ to parts in 10^ from the momentum-angle of the decay positron. The 



polarization of the muon from pion decay is P^. The SM predicts that p = 
and ?? = 0; deviations from these predictions would indicate new physics. 



3/4, = e = 1, 



2. Experiment 

The experiment used the M13 beam line at TRIUMF in Vancouver, Canada. Positive pions 
decaying at rest at the surface of a carbon production target produced highly polarized positive 
muons with a momentum of 29.792 MeV/c. The M13 channel momentum selection was set at 
29.6 MeV/c with a resolution of 0.7% to select these muons and guide them to the spectrometer. 
A thin metal foil acted as a stopping target for the muons; this was placed at the center of a 
symmetric array of wire chambers within the bore of a solenoid that produced a highly uniform 
2.0 T magnetic field known to three parts in 10^ (see Fig. [T]). The particle identification relied 
mostly on three modules of four proportional chambers (PCs) each. PC modules were installed 
at each end of the spectrometer. The third module used the stopping target as a cathode foil 
in the center of the four PCs and was installed in the center of the detector stack. Muons that 
stop in the target were selected using the PCs of this target module. The decay positron helices 
were tracked by 44 drift chambers, and their trajectories were later reconstructed to determine 
the positron's initial momentum and angle. The wire positions were known to five parts in 10^, 
providing a high reconstruction resolution of 60 keV at a positron energy of 52 MeV. The wire 
chambers were low mass to reduce multiple scattering and to allow the muons to reach the target 
since it takes only about 1 mm of water equivalent to stop muons at 29 MeV/c. Further detail 
on the apparatus can be found elsewhere [2]. 
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Figure 2. Experimental momentum-angle spectrum. The detector response is included in this 
spectrum. 

The muon decay parameters were measured by comparing the positron momentum-angle 
spectra (see Fig. [2|) from the data to a GEANT3.21 simulation that was subjected to the same 
analysis. In this way the detector response and reconstruction biases were accounted for within 
the simulation. Hidden values of p, 6 and ^ were used in the simulation, and these were not 
revealed until the corrections and systematic uncertainties had been evaluated on the difference 
in decay parameters between the data and simulation spectra; this technique provided a blind 
analysis by exploiting the spectrum's linearity in p, and PJ^^S (see Eqs. (I3]),(jl|)). 

Special data validated the positron physics in the simulation using muons stopped close to the 
entrance of the detector. In this configuration the decay positrons traversed the whole detector. 
The corresponding tracks were independently reconstructed in each half of the detector, before 
and after crossing the stopping target. The reconstruction efficiency was measured from this 
special data by counting the number of tracks reconstructed in one half of the detector but not 
in the other. 

3. Improvements for the final measurement 

An initial and an intermediate measurement of the decay parameters were already performed by 
the TWIST collaboration on data taken in 2002 [3l|4] and 2004 OE]. Final data were acquired 
in 2006 and 2007, with a higher quality muon beam and a threefold increase in statistics. 

A new technique was developed to measure the space-time relationships (STRs) used to 
convert the drift times into drift distances in the drift chambers [7]. These improved STRs 
also corrected reconstruction biases. Each drift chamber was calibrated independently; this 
accounted for small differences in construction and response. The beam line was upgraded 
to correct an undesirable muon beam vertical deflection of about 1.0 cm. The beam was 
steered onto the symmetry axis of the solenoid, which reduced the uncertainty in simulating 
the depolarization of the muon. The long term stability of the beam was monitored using its 
average position measured by the wire chambers. Muons were stopped in both an Al and Ag 
target foil (previously only an Al foil was used) in two separate datasets. This allowed for the 
study of systematics depending on the properties of the target material. 



Table 1. Systematic uncertainties for each decay parameter in units of 10~^ from the bhnd 
analysis. The depolarization and background muons systematic uncertainties affect only the 
muon polarization P^. 
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Beam stability 
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Radiative corrections 
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Total 


2.8 


2.9 


+16.5 
-6.2 



4. Systematics uncertainties 

The systematic uncertainties from the blind analysis of the final measurement are summarized 
for each decay parameter in the Table [H 

Most systematic uncertainties are evaluated by altering the component source of uncertainty 
in the simulation or in the analysis. The momentum-angle spectrum created by this modified 
simulation or analysis is then fitted against the unaltered spectrum. The alteration is typically 
many times greater than the measured uncertainty on the modified component in order to 
increase the sensitivity. For this reason the difference in decay parameters is scaled down 
to match the uncertainty on the modified component. This scaled difference is used as the 
systematic uncertainty on the decay parameters. For example the bremsstrahlung production 
rate uncertainty (the dominant uncertainty in the group "positron interactions") was evaluated 
by generating a simulation with the bremsstrahlung production rate multiplied by a factor of 3. 
The bremsstrahlung rate was measured in the standard simulation and the experimental data 
using the topology of the events to identify events with a bremsstrahlung being emitted. The 
simulation and the data bremsstrahlung rates differ by 2.4%. Therefore the difference in decay 
parameters between the altered and the unaltered simulations was multiplied by (1.024-l)/(3-l) 
to provide the corresponding systematic uncertainty on the decay parameters. 

The dominant systematic uncertainty for the parameter comes from the depolarization 
undergone by the muons as they enter the 2.0 T tracking magnetic field. A mismatch in 
the depolarization in the fringe field region between the simulation and the data leads to a 
mismatch in the muon polarization at the time of decay. This creates a systematic bias in 
the determination of and in the measurement of P^S,- The accuracy of the simulation of 
the depolarization was evaluated by modifying the position or size of the experimental muon 
beam and verifying that the simulation could reproduce the resulting change in polarization. 
This evaluation indicated that the simulation underestimates the depolarization which leads to 
an asymmetric uncertainty. 
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Figure 3. Blind analysis results plotted along with the previous TWIST and pre- TWIST 
measurements. 



5. Results 

The TWIST collaboration agreed on the list and values of the systematic uncertainties and 
corrections before revealing the hidden parameters of the blind analysis. 
The results of the blind analysis are 

p = 0.74991 ± 0.00009 (stat) ± 0.00028 (sys), (5) 
S = 0.75072 ±0.00016 (stat) ±0.00029 (sys), (6) 
P;^ = 1.00083 ± 0.00035 (stat)l°:°°J^^ (sys). (7) 

The parameters p, 5 and are respectively 0.3, 2.2 and 1.2 standard deviations away 
from the predictions of the standard model. All the parameters are consistent with the previous 
measurements from TWIST and from experiments prior to TWIST [SlIHKlO] (see Fig. [3]). 

The spectrum asymmetry Aep at the positron kinematic end point is given by 

Aep = therefore < 1. (8) 

P P 

However the results from the blind analysis give: 

5^ - 1 = (1921^^^) X 10"^ (9) 

which corresponds to 2.9 standard deviations above the physical limit of one for the four-fermion 
interaction model. At the present time we assume that there is a systematic uncertainty or 
correction that we haven't identified. For this reason the blind analysis results are not considered 
final and could be subject to change. 

6. Theoretical implications 

The blind analysis results can be used to put stringent constraints on new physics. It is important 
to emphasize that these constraints like the blind analysis results given above are not to be 
considered final. 

In Left-Right Symmetric (LRS) models the right-handed current is suppressed but not zero. 
An additional heavy right-handed W-boson (Wr) is introduced to restore parity conservation 



at high energies pjj. In these models, the left- and right-handed gauge boson fields are given 
by: 



Wl 

Wr 



Wl cos C + W2 sin C 
e^'^(-TyisinC + W2C0SC) 



(10) 
(11) 



where u is a. CP violating phase. The TWIST result for p allow for model-independent 
constraints on the mixing angle (C) between the Wl and Wr and on the mass 771-2 of the W2 
mass eigenstate. No assumptions on the left and right couplings, CKM matrices, or on the 
CP violation are made. The pre- TWIST limits from muon decay were ((7r/(7l)|C| < 0.06 and 
{or/ 9l)'i^2 > 400GeV/c^. Our preliminary results improve these limits to {gR/gL)\C\ < 0.02 
and {gR/gL)m2 > GSOGeV/c^. 

7. Conclusions 

The blind analysis of the final data from the TWIST experiment reached the precision goal of 
a few 10~^ on the measurement of the decay parameters p, 6 and P^^. The product Pp^6/p is 
2.9 standard deviations above the physical limit of one defined by the four-fermion interaction 
model used. The present results are therefore not final and the possibility of a missing systematic 
uncertainty or correction is being investigated. 
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